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Abstract 

We propose the method of statistical description of broad class of dynamic 
systems (DS) whose equations of motion are determined by two state depending 
functions: 1) "energy" - the quantity which conserves in time and 2) "entropy" 
- the quantity which does not decrease in time. It is demonstrated that the 
behavior of such systems in the equilibrium state reduces to the thermodynamic 
lows in particular the Le Chatelier principle is satisfied and so on. Taking into 
account the interaction system of interest with ergometer - the device which 
continuously measures its energy one can possible to find the system distribution 
function in arbitrary non-equilibrium stationary state (NESS). Some general 
relations for mean values of certain quantities in NESS which can be compared 
with experimental data are obtained. 
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1 Introduction 



The dynamic equations method is the base mathematical tool for investigation 
of the behavior of complex systems. Such method of description is widely used 
not only in classical physics where it appeared for the first time but also in 
chemistry, biology, ecology and other sciences. The types of equations describing 
sufficiently broad classes of systems (and the processes taking place in them) 
on the one hand and having some general features which do not depend on the 
nature of the studying objects on the other hand are especially important. In 
this connection we should point out such classes of equations:!) Hamiltonian 
ones and 2) equations of gradient type. Hamiltonian equations of motion have 
the form: ^ = |^;^ = —§^, where H{xi,p.j) is the Hamilton function 
of the system. They are used for describing the evolution of conservative DS 
with even-dimensional phase space whose behavior is reversible in time. For 
describing the nonreversible systems with the strong dissipation the equations 
of gradient type are usually used: ^ — where '^{xi) - Lyapunov function 

of the system, the quantity which monotonically decreases in time. 

Let us now call attention to the study of evolution of closed system with the 
large number of degrees of freedom. It is well known that such system comes to 
the equilibrium state without fail at the condition i — > cxd (zero law of thermody- 
namics). In his seminal paper of 1865, R. Clausius stated the laws determining 
the behavior of such systems in the form of the following two propositions (see 
[I],p. 194): 

1) Die Energie der Welt ist constant (1) 

2) Die Entropie der Welt strebt einem Maximum zu 

We'll consider the class of DS which satisfy to this Clausius conditions. In 
this connection we want to ask two questions: 1) What equations of motion 
with necessity provide such behavior of the system? 2) What results related to 
the statistical behavior of such systems in and out of equilibrium can be ob- 
tained from such equations of motion? The attempt to answer these questions 
is the main goal of our paper. Note that in phenomenological approach which 
is used in this paper the terms "energy" and "entropy" will be used only as 
conventional labels for two quantities satisfying the Clausius conditions. Now 
we describe briefly the content of the paper. It consists of five sections. In Sec.l 
we introduce the notion of quasithermodynamic system (QTS) and give the con- 
crete recipes to construct the equations of motion for such systems. In Sec. 2 on 
the base of these equations we investigate the conditions of equilibrium of two 
QTS. We prove that in the case when an interaction between two QTS is ab- 
sent our conditions of equilibrium are reduced to the standard thermodynamic 
relation: = Besides using the condition of equilibrium state stability 
the Le Chatelier principle for QTS is derived directly from equations of motion. 
In Sec. 3 we propose the method of statistical description of QTS behavior. In 
order to avoid the appearance of singular distribution functions we introduce 
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the interaction between the system of interest and the ergometer- device which 
continuously measures its energy. Such method gives us the possibihty to for- 
mulate the simple Fokker-Planck equation for QTS distribution function whose 
stationary solution can be obtained. This solution can be expressed via state 
functions S{xi) and E{xi) and additional parameter /3. In the limiting case 
of equilibrium state the solution obtained reduces to the ordinary Gibbs dis- 
tribution: f{xi)~ exp{—PE{xi)). In Sec. 4 we investigate the non-equilibrium 
stationary states of QTS more detail and formulate some general relations for 
them. In the last part of the paper we summarize the results and briefly discuss 
the possible applications of them. Now let us pass to the detailed account. 



2 Section 1 

We define as quasitermodynamic (of N order) the DS whose state is described 
by N variables Xi^..,Xn and whose behavior in time satisfies two Clausius con- 
ditions (1). Meanwhile we suppose that variables Xi^..,Xn unlike from coor- 
dinates and momenta define the "coarse -grained" description of macroscopic 
system evolution and therefore the number of them can be arbitrary. Let us 
establish now the form of equations of motion for given state functions: en- 
ergy H{xi^..,xn) and the entropy S{xi^..,Xn) providing the realization of the 
distinctive properties of QTS : 

dff _ dff ^ _ ^ 

dt dxi * ' 

dS dS ■ 

(In relations (2), (3) and later on we use the usual agreement about the 
summation on repeating indexes). 

We start from the simplest case of QTS of second order (A'' = 2) and sup- 
pose that required system of equations of motion is an autonomous system 
of differential equations: ^ = ^i{xi,X2)- By direct verification one can be 
convinced that taking functions ^i{xi,X2) in the form: $j = eik§^ {S, H}, 

where Sik = ^ ^ ^ is antisymmetric tensor of 2 rank and {S, H} is 

Poisson bracket of functions S{xi,X2) and H(xi,X2), i.e. by the definition 
{S,H} = - we obtain the equations of motion 

from which one can easily obtain properties (2) and (3). 
Similarly in the case TV = 3 one can be convinced that the system of equa- 
tions of motion with required properties can be presented in the form: 

dXi dH 
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where Siki is the antisymmetric tensor of 3 rank and 



dS dH 

Ai=ei^^^ — . (6) 
In addition we give the form of equations of motion for the case = 4: 

dxi dH ^ , 

at OXk 

where Sikim is the antisymmetric tensor of 4 rank and 

Aim = eimnp-K—-K—- {°) 

\y JU Y}^ *Xy ip 

The generahzation of equations of motion in the case of arbitrary N is ob- 
vious. We wiU not discuss here the rigorous mathematical justification of the 
uniqueness of equations obtained. About this problem one can see the paper [2] 
and references in it. In this paper we are interested in the properties of QTS 
implied from equations of motion (4), (5), (7) and the possibility of statistical 
description of such systems. 



3 Section 2 

We start our analysis of equations of motion for QTS from the simplest case 
^ = 2. 

Let us rewrite the equations (2) into another, more convenient form. For 
this purpose we introduce the "phase' ip which is variable canonically conjugated 
to the energy i.e. {(p,H} = 1. Now one can easily see that the system (2) is 
equivalent to the following system of equations for H and ip: 

dt ' dt d^' ^ ' 

where QTS entropy S now must be expressed via ip and H. It follows from 
(9) that in the case when S{ip, H) has maximum on ip there is the line i.e. the 
one-parameter set of equilibrium points ip = Pq{E) corresponding to the solution 
of equation ^^'^q^^^ — when ip — ipQ^E). Similar equilibrium condition is valid 
in general case also. It is easy to show that equilibrium conditions for QTS of 
A'^-order reduce to the next system of equations (it follows from equality to zero 
the right parts of (4), (5), (7)): 



dS 


as 


as 


dxi 


dx2 


axM 


IK 


~ dH_ ~ ■ 


■■ an 


dxi 


8X2 


axM 



A. (10) 

Now let us consider the composite system consisting of two subsystems: A 
with variables xi^..^xni and B with coordinates yi,..,yN2- Besides we suppose 
that subsystems A and B are non-interacting i.e. relations H = Ea{xi^.., xni) + 
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EB{yi.--, yN2) and S = Sa{xit.,Xni) + SB{yi,--,yN2)- In this case equilibrium 
conditions reduce to the form: 



OSa OSa OSb dSE 

(11) 



dxi 


dxNl 


dyNi 


dxN2 


OHa ~ ■ 


■■ dHA 




■ OHb 


dxi 


dxMl 


dvNi 


dxN2 



Equations (11) can be satisfied by setting Sa = Sa{Ea) and Sb = Sb{Eb) 
after that equations reduce to standard condition of thermodynamic equilibrium 
of two subsystems: 

dSA^dSB^ , . 

Now wc demonstrate how using only cquihbrium equations (10) one can 
deduce the Le Chatelier principle. We remind that Le Chatelier principle claims 
that response of the system in the equilibrium state on the external influence 
always results in the decreasing of this influence. For proving this principle 
in the simplest case of 2-order QTS with variables x and y let us consider 
three states of such system: I - {xo,yo) - the initial equilibrium state, II - 
(xq, 2/0 + Sy) - nonequilibrium state enforced by small influence on the variable 
y, III - {xo + Sx, yo + Sy) - the final equilibrium state. It is obvious that the final 
equilibrium state III differs from initial equilibrium state I. Because the energy 
of the system have been changed. The first condition wc write down is equality 
of energies of states 11 and 111: Hu — Hm — H{xo, yo) + SH, where 

Besides both states 1 and 111 are equilibrium states the following conditions 
must be satisfied: 



(ii)o ^ {^y)o 

(m fdH\ 

{Ox)o [g^)^ 



(14) 



and 



"I" SxxSx + SxySy Sy + SyxSx + SyySy 



(15) 

Hx + HxxSx + HxySy Hy + Hy^Sx + HyySy 
where for simplicity we used notations: Sx = (ff)Q = (if)(xoyo)' '^^ ^ 
and so on. We can write Sx = ASy and Sy = BSy, where A and B are 

yet unknown functions of Xo and yo. Equations (13) and (15) now are reduced 
to the form: 

Hy = HyB + HxA, (16) 
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{Sxx-^ + SxyB) Hx Sx {Hxx-^ + HxyB) 

Hi " 

(^Syx-A. + SyyB) Hy Sy (^HyxA + Hyy B) 

~ IP ■ 

V 

Wc can simplify the form of equations (16) and (17) by using the quantity 
F = S — XH. Then taking into account (14) for unknown functions A and B 
we obtain the following system 



H^A + HyB = Hy, f^^^)A+(^-^)B = 0. (18) 

Solving the system (18) we obtain 

H H 

~ {HxFyy HyFxy^ 7 B — ^ {^HyF^X -^xFxy) 5 (-^^) 

where A = H^Exx - 2HyHxFxy + H^Fyy. 

Note that equilibrium stability conditions for state I require inequalities 
F^x < 0, Fyy < and F^^Fyy — Pxy > ^ ^'"^ valid. In particular this case follows 
that A < 0. 

Now for proving the Le Chatelier principle it is necessary to introduce the 

measure of QTS reaction on the exterior influence. As such measure it is natural 

to choose the magnitude of thermodynamic force Y = — [ ^ J that is equal 



dy ^ 

to zero in unperturbed state I. Then for state II we have following relation: 
AF = Yii -Yi = {^%)jy = -FyySy and for equilibrium state III AY = 

Yni -Yj = (||)„ ^ + (^)^ ^ = _ ^FxyA + FyyB) 6y. 

The Le Chatelier principle will be valid under the condition 

FxyA + FyyB ^ Fyy. (20) 

Substituting A and B from (19) into (20) after some simple algebra we came 
to the inequality FyyH^ + HyF^y ^ 2HxHyFxyFyy, which is obviously valid. 
Thus the validity of the Le Chatelier principle is proved. 



4 Section 3 

Now let us pass to the statistical description of QTS and for simplicity we start 
from the system of 2-order. Our main goal is to determine the distribution func- 
tion for nonequilibrium states of such system. Note that under pure dynamical 
description the distribution hmction f{xi,t) must satisfy the continuity equa- 
tion: ^ + div{fxi) — 0. Because of attractors existence in QTS such function 
necessarily must be singular. To get the smooth distribution function we will 
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use the following method. We introduce the interaction of QTS of interest with 
ergometer - large exterior system which measures its energy From the definition 
of ergometer it follows that during the measurement process mean energy of 
QTS system remains the same. As we have showed in |3] such method results 
in Fokker-Planck equation for f(xi,t) of the form 

where diffusion tensor Dik representing the influence of ergometer on the 
system behavior must be chosen in the form: 

dH dH 

^ik = £il£kmT: — , (22) 

dxi dx,n 

Here H{xi, X2) is the energy of QTS. It is obviously that tensor Dik defined 
in such a way possesses the symmetry property: Dik = Dki and the nonnegative 
property: Du > 0, D22 > and D11D22 — Df2 > 0. Having the Fokkcr-Planck 
equation (15) in hand we can search its stationary solutions corresponding to 
equations of motion: Xi = eik-§^ {S, H}. It is easy to show that all such 
solutions can be represented in the form: 

/(a;i,:E2) =A(ff)e^(-^'-^), (23) 

where S is the entropy of QTS and A{H) is arbitrary function of its energy. 
We define function A{H) from the following physical reasons. First of all let us 
pay attention to the fact that QTS entropy is defined up to arbitrary function 
of energy because the transformation: S'(a::i,a;2) S{xi,X2) + C{H) does not 
change the form of equations of motion. Besides we must require the distribution 
function (17) transforms into well-known Gibbs distribution in the limiting case 
of equilibrium state. It is easy to verify these conditions impose by the only 
possible choice of distribution function in the form: 

f{xi,X2) = Are^(^i'^^)-'^='(^?'^^)-/3^^(^i,^2)^ (24) 
where TV is normalizing coefficient defined by the condition: 



J f{xi,X2)dxidx2 = 1, 



and Seq{xi,X2) is the entropy in the equilibrium state what is the function 
of the QTS energy and /3 is the parameter having sense of inverse temperature 
of system in the equilibrium state. The relation (24) for distribution function of 
QTS is valid both in equilibrium and nonequilibrium states and it is a main re- 
sult of this section of the paper. Note that relation of type (24) for distribution 
function can be derived also on the base of information-theoretic Jaynes method 
[4] but we suppose these two methods are essentially different. The main differ- 
ence that in our method the " entropy" is actually the dynamic quantity defined 
by equations of motion. So it can be used for describing both equilibrium and 
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nonequilibrium states of QTS. Also it is worth to note that for deriving the 
distribution function (24) from Fokker-Planck equation the Markov approach of 
system evolution is required. This fact is ignored in Jaynes method. 



5 Section 4 

In this section using the explicit form of distribution function (24) we find 

the number of relations which are connecting the different characteristics of 
QTS. First of all we point to simple relation between relaxation time t{E) and 
entropy of QTS in the equilibrium state. For this purpose we write down the 
equation of motion for phase close to equilibrium line (p = (Pq{E). In the lowest 
approximation it has a form: (p = — ■ Comparing this equation with 

(9) we derive: 

t{E) = --. ^ . (25) 

It is obviously not far from equilibrium line QTS entropy may be written in 
the following form up to terms of second order: 

5(^,^) = 5e,(i;)-fc^p!. (26) 

Comparing (20) with (18) we derive the distribution function for NESS close 
to equilibrium: 

f{<p,E) = Ne (27) 

We denote by Aip the deviation of system phase from its equilibrium value: 
Aip = if — ipq{E). Let us evaluate the mean value ^ ^t{e) ) NESS. From the 
definition of mean value we obtain 



/ (A^)^ \ _ / ^e-m-^-d^dE _ J ^)e-0-dE _ ^ ^^^^ 



\t{E)/ j^-(-^y-0E^^^^ J,/^e-f^^dE 
Relation ( ) = 1 is valid close to equilibrium state as one can see 

\ ■^(^) / NESS 

from its derivation. Now we derive the relation which is valid not only close 
to equilibrium. Let us consider the mean value (^ ^^^qe^"^ ) • Using the 
definition of mean value and integrating by parts in numerator we obtain: 

/d{S-S^\ _ J 'J^e^-s..-P^d^dE _ 

\ dE / je^-^'^-P^difidE ^ ' 
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The obtained relation (29) connects mean value of '^^'^g^'"''* in nonequilib- 
rium stationary state with inverse temperature of equilibrium state of the same 
mean energy. And finally let us set the connection between the relaxation time 
t{E) and inverse temperature /3. For this purpose we consider the value 

oo 

Z{P) = j y/^e-'^^'dE. (30) 



As follows from (22) the function ^T{E)e~^^ has the sense of non-normalized 
distribution function of QTS close to equilibrium. Now we find the mean energy 
of the system - U. By the definition we have: 



We suppose t{E) to be monotonically increasing function of E and also 
suppose that the function ^jT{E)e~^^ possesses the sharp maximum on E. 

Then we can find the value Z{P) using the saddle point method: Z{(i) w 
T{U)e~^^ , where mean energy U is defined by the condition of integrand func- 

ax ' 

tion maximum in (30) — j3\/T{U) = 0. From that we find the required 

connection between r and /?: 

/3(t^) = ^^lnr(C/). (32) 

Note though we investigated only the simplest case of QTS of 2-order nev- 
ertheless all results obtained in Sec. 3 and Sec. 4 may be generalized with the 
small modifications to the case of arbitrary A^. For example for A^ = 3 the dif- 
fusion tensor corresponding to the interaction with ergometer must be chosen 
in the following form: Dik — Siim£kin-§^§^- One can see the form of station- 
ary distribution function (24) derived as the solution of Fokker-Planck equation 
(21) remains the same. We give also the analog of (32) for QTS of A'-order 
without derivation. In this case system has the collection of relaxation times: 
Ti{E),T2{E), ...tn-i{E). They are connected with the inverse temperature of 
the equilibrium state by the relation: 



JV-l 



In n r^{E) 



(33) 



6 Conclusion 

Let us sum up the results of the paper. Our main goal was to demonstrate how 
using Clausius conditions (1) one can successively construct classes of DS of 
increasing complexity exactly satisfying to these conditions. We demonstrated 
also that properties of such systems in the equilibrium coincide with properties 



9 



of usual thermodynamic systems. That's why the term " quasithermodynamic 
system" was proposed. Besides using the interaction of QTS of interest with 
ergometer we could give the statistical description of nonequilibrium stationary 
states and obtain universal distribution function for them. Note though we 
investigated the case of classical systems only the proposed method can be 
extended to the quantum realm. Speaking more accurately the equations of 
motion for QTS (4) admit the semiclassical quantization, for example by using 
the method of paper [5] . Detailed investigation of quantum- mechanical variant 
of QTS description will be given in future publications. Concerning to the 
possible application of obtained results we must point out that they are highly 
diverse and may be applicable not only in physics. We bring here only one such 
example concerning to the chemical kinetics. It is known ( see e.g. [6]) that for 
closed chemical systems where the chemical reactions not higher then second- 
order on concentrations of reagents Xi ^ take place one can write down the 
following system of kinetic equations: 

= /3,feXfe -f JikiXkXi. (34) 

We suppose the absence of autocatalytic reactions in such system and con- 
dition p^Xi — M where is molecular mass of reagent Xi is satisfied. Then 

i 

it is easy to demonstrate that all such systems are quasithermodynamic. Let 
us illustrate this statement on the simplest case of two reagents xi and X2 with 
coinciding molecular masses: Pi = P2- 

In this case kinetic equations have the form: 

= -axi + (3x2 ~ JX1X2, = axi- 13x2 +1x1X2- (35) 

One can easy see by the direct verification that using functions 

H = x,+X2, S^^- ^ —. (36) 

we can rewrite the system (35) in a standard form: = Sik^^ {S^H} . 
This fact means that all results of the paper, in particular the possibility of 
description of reagents concentrations fluctuations, become applicable to the 
closed chemical systems. From this example we can see also that the nature of 
functions H and S is defined exclusively by character of the considered system 
that is by its equations of motion and may be essentially different from their 
thermodynamic analogs. 



Author would acknowledge L.A. Pastur for the discussion of the results of 
the paper and valuable comments. 
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